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ABSTRACT 

Advances in digital computer technology have led to employment of 
the computer as a control device. Although many theories for utilization 
of the digital computer as the compensator in a position feedback control 
system have been proposed, there is little documentation of digital 
control performance in a hardware environment. This thesis is an 
investigation of the performance of a real system under sampled-data 
control. The variable gain amplifier sampled-data control theory is 
first simulated and then tested on a real system. The hardware and 
software developments necessary for implementing this control theory 
are discussed in detail. Additionally, this method is incorporated in 
a time-sharing computer program for controlling many systems simultane- 
ously from a single computer. Results of the tests are presented and 
evaluated. 

The authors wish to express their appreciation to Professor Robert 
D. Strum of the U. S. Naval Postgraduate School for his guidance and 


assistance during this investigation. 
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CHAPTER I 


iL SintroducrLen 

During the past 15 years, advances in digital computer technology 
have developed the digital computer into an invaluable tool for all 
branches of the scientific world. In almost every scientific discipline, 
the capabilities of the digital computer for vast computational tasks, 
data processing, and simulation have been profitably employed. In 
addition to these uses, some engineers in the automatic control field 
have considered the digital computer, operating in a real time mode, as 
a tool for improving control techniques and capabilities. It is in this 
area of computer application that the authors are extremely interested. 

There are many theoretical expositions on sampled data control 
techniques, and some include results from the digital control of 
simulated plants. However, with the exception of the technical reports 
of the Navy Electronics Laboratory EE | the authors were unable to 
find any information on sampled-data control of an actual plant. It is 
the opinion of the authors that unless the theory is tested in a hardware 
environment, the results are academic because of the difficulty in 
accurately simulating the innenene non-linearities in a real system. It 
was decided at the outset, therefore, that a prime goal of the thesis 
was to test the theories of sampled-data control on a real plant. 

A second area in which data was lacking involved one of the primary 
advantages claimed for the sampled-data control techniques: the cap- 
ability to control many systems simultaneously from the same computer. 


Although implementing such a multi-~control system is primarily a computer 





programming task, certain basic techniques must be developed to insure 
flexibility and adequate performance. It was felt that the lack of 
documentation in this area was a serious shortcoming. Thus, a second 
goal of this thesis was to develop and implement techniques for multi- 
system sampled-data control. 

Having established the basic goals of the thesis, the authors made 
a decision which essentially limited the scope of the thesis. Throughout 
the literature, there are many proposed theories for sampled-data control 
systems. It was decided to select only one of these theories as the 
basis for the laboratory experimentation. Rather than testing many 
theories in a less rigorous fashion, the authors preferred to make 
exhaustive hardware tests on the basis of one theoretical solution to 
the control problem. In general, the difficulties encountered in one 
theoretical approach would be comparable to those of another, and the 
techniqués developed may be extended to other theories. 

After examining many of the proposed theories, the authors, in 
conjunction with their advisor, Professor Robert Strum, selected the 
theory proposed by Benjamin C. Kuo. The basis of this theory is the 
concept of a variable gain amplifier in the control loop. The theory 
provides for minimum time, ripple free response to deterministic inputs. 
The prime reasons for selection of this theory as a basis for the 
laboratory experimentation are the theoretical advantages it offers, 
and the relatively simple tasks the digital computer is required to 
perform. 


In considering the laboratory portion of the thesis, the authors 





had one more basic decision to make. In selecting a digital computer 
for the experiments, the facilities at the Postgraduate School offered 
a choice between a general purpose computer, the Control Data Corporation 
1604, and a small data processing computer with an auxiliary arithmetic 
capability, the CDC-160 in conjunction with the CDC-168 arithmetic unit. 
Due to the fact that any practical control method must not only perform 
well, but also must be economically feasible, the smaller, less expensive 
160 computer was selected. It is conceded that the larger computer could 
more efficiently perform the required controlling tasks, but the authors 
preferred the more realistic practical approach. This led to a secondary 
goal of the thesis, to be able to comment on the general feasibility of 
sampled-data control systems on the basis of the results obtained. 
1.2 Digital Controller Design 

Before proceeding with work on the real plant, the authors decided 
to select a technique for the design of digital controllers for sampled- 
data control systems. A design theory that would give accurate results 
and would be readily adaptable to the existing hardware was desired. It 
was determined that the digital controller design described in Section 
9-9 of Analysis and Synthesis of Sampled-Data Control Systems by 
Benjamin C. Kuo would be well suited to this problem. [ 3 ] Kuo's design 
theory makes use of state space and state transition techniques. It 
also enables derivation of a D(Z) that will yield a minimum time, ‘ripple 
free response to the designated deterministic input (i.e., "deadbeat 
response''). This approach to digital controller design describes the 
digital controller as a variable gain amplifier which outputs certain 


values each sampling period in order to control the continuous system 








in the desired manner. It is easily adaptable to implementation on a 
digital computer, and the use of state variables makes the manipulation 
of the mathematics involved relatively easy. 

To explain this "variable gain amplifier'’ design theory it is best 
to look at a sample problem. Consider the sampled-data control system 


of Figure 1. The variables xy and X, are the state variables for the 


controlled process where x, = x h(t) is the output of the zero-order 


2? 
hold and e(t) is the actuating signal. Note: e(t) = r(t)-c(t)5x,=c(t) 


i} 
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X | Zi Oot h¢t 1 
sis s(st+1) x,=c(t) 


Le G(S) 
Figure l 

To determine a D(Z) for "deadbeat response" to a given input: 

(1) Draw the state transition flow graph of the system with the 
digital controller represented by a variable gain amplifier 
Koch ey Ken) — hk.) o/ etn 

(2) For a "deadbeat response" the system error must be zero for 
t >nT where n is the smallest possible integer. 

i.e., xy (nT) = r(nT) 
Xo (nT) = x, (nT) ere ee = 0. unit step 
Xo (nT) = 1; eer) aio EE) = 0. unit ramp 
From this information the variable gains (K. 's) and the minimum 


mM Can be determined. 





(3) Compute the h(0')., h(T'), .. h(t) using the Ky Ne 
determined and recalling that e(kT")=r(kT)-x, (kT). 
(4) Calculate D(Z) using the relationship D(Z)=H(Z)/E(Z). 


Returning to the example, let us draw the state transition flow graph 


of the system. 





Figure 2 
Using Figure 2 and applying flow graph techniques and the inverse 


Laplace transformation we arrive at the following: 


xy [ cctay | 
Xo [at | 


Now let us assume a unit step input and a sampling period of one tenth 


[1-K, (t-14e77) ] x, (KT)+(1-e7")x, (kT)+K, (T-1+e™") x (kT) 


-K,(1-e")x, (kT) +e" "x, (KT)+K, (1-e” ")x (KT) 


of a second. ( r(kT) = 1.0, T=0.1) 


xy | cet) | 
Xo | (k+1)T ] 


For k=0 and x, (0) = x, (0) fo. 


(1 - 0.005K, )x, (kT) One o> x» (kT) te 0.005K_ 


-0.095K_ x, (kT)+ 0.905x, (kT) + 0.095K,. 





x, (T) 0.005K. 


x,(T) = 0.095K 


From the constraints placed on the problem (i.e., ee = r(nT) and 
x, (nT)=0. ) it may be seen that the solution can not be obtained from 
the above. In practice this means that the output of the controlled 
process cannot be made equal to a stép input in one sampling period. 


We must therefore proceed to the second sampling period and let k=1. 


Doing this and setting x De 1.0 and x (2T) = 0. yields: 
x, (2T) = 0.014K - 0.000025K K, + 0.005K, £ 1.0 
l O fem! i 
x,(2T) = 0.086K_ - 0.000475K K. + 0.095K., 2 0.0 
2 O cal 1 
from which 
K = 105.6 ; K, = -203.0 
O il 
therefore 
“ + 
e(0 ) = r(0)-x, (O)=1 ; e(T =r (T)-x, (T) = 0.472 
OOD), & ese - h(t’) = Ke (I) = -95,8 
and 
D(z) = n(Ot) + h(tt)z"* _ 105.6 - 95.8 z+ 
AGE ye it OA 


The above example demonstrated the technique of designing a digital 
controller for a sampled-data control system using the variable gain 
amplifier concept. It should be observed that in implementing this 
controller using a digital computer, the D(Z) does not have to be actually 


calculated. In fact, only the Ks need be entered in the computer. 





This proves to be quite advantageous since for a given sampling period 
the K's are only a function of the continuous system gain and the G(S). 
Thus for a digital computer controlling a number of plants the Ks 
could be pre-computed and stored in the computer ready for use. These 
variable gains when multiplied by the sensed error and output at the 
sampling rate will provide the desired system response to the designated 
input. In fact, the digital controller designed above will yield the 
desired "deadbeat response" to any step input. It must also be pointed 
out that the design given is only for a step input and will not give 
optimum response for a ramp or parabolic input. 

The design problem as stated and solved above was for the control 
of a type one system in response to a step input. [In studying this 
design technique a digital controller was also designed for the same 


type one system (G(S)= ) for "deadbeat response" to a unit ramp 


1 
s(st1) 
input. With respect to a ramp input, optimum response was obtained in 
three sampling periods. It was also noted that once zero error had 
been achieved a constant output of 1.0 was needed from the controller. 
This was required because a type one system has an inherent Steady state 
Paeor in response tO a ramnpeinpuc. In this context it is important te 
realize that with the use of a digital controller a continuous system 
with an inherent steady state error can be made to respond with "deadbeat 
response" and no steady state error to any particular deterministic 
input. (See Appendix I for details.) 

In Appendix I a digital controller is also designed for a type 


zero system (G(S) = eS EOE) ) which has an inherent steady state 





error to both step and ramp inputs. With the designed controller, these 
errors were reduced to zero. 


With regard to type two systems some difficulty was encountered. 


For a system of the type G(S) = = a digital controller was readily 
S 
designed; however, for a system of the type G(S) eee the mathe- 
s (stl) 


matics involved in evaluating the K,'§ proved to be cumbersome. An 
iterative solution to the non-linear simultaneous equations was used 

in evaluating the KS: This proved to be feasible. Further investiga- 
tion into the problem was not carried out since the authors wished to 


direct their efforts to the control of the real system. 


1.3 Simulation and Results 

After completing the design of a digital controller for a 'dead- 
beat response" to a unit step input, it was decided to simulate the 
system response. This was done to provide a check on the computations 
and to verify the design theory itself. The simulation programs used 
were written in Fortran 60 and employ a library routine known as INTEG l. 
INTEG 1, written by Dr. J. R. Ward, provides a fourth order Runge-Kutta 
solution to ordinary differential equations and was well suited to the 
simulation problem encountered. (Appendix II contains a complete 
simulation program. ) 

The graphical results to the simulation are on the following pages 
and do verify the design theory as well as the computations involved. 
(Appendix I contains graphical results for a unit ramp input as well 
as results for a type zero system in response to both unit step and ramp 


inputs.) 





It should be mentioned that a certain degree of accuracy is needed 
in the determination of the h's used. In some simulations, the rounding 
of the h value to two decimal places resulted in a response that was 


slightly less than optimum. 
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CHAPTER [II 
2.0 Investigation 
The investigation conducted during the course of this thesis was 
developmental in nature. Using the variable gain amplifier control 
theory discussed in Section 1.2, techniques for accomplishing the 
sampled-data control of a real plant were developed. The investigation 
also included the development of a computer time sharing system through 


which many systems can be simultaneously controlled. 


2.1 Theoretical Design 





Figure 6 


Figure 6 shows the general block diagram of the sampled-data 
control system that was used. The system breaks down into two main 
components which are the continuous plant and the D(Z), sampler, and 
zero-order hold. The continuous plant is described in Appendix III. 
The CDC-160, CDC-168, and A/D, D/A converters function as the sampler, 
D(Z), and zero-order hold. This operation will be treated in detail 
Hi OOCE LOM tee = ly. 

The digital controller, D(Z), was designed to give a "deadbeat 


response" to a step input using the theory described in Section 1.2. 
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It was decided to solve for the K's in general terms initially and 


then substitute the known G(S) determined in Appendix III. 


ae ee le 
S(sta). 


Figure 7 is shown for a generalized G(S) = 





Figure / 
From the flow graph we have 
i KA 1 KA 
me ()-[2 - 2} (t )t+ — )SO=xA(t_)+ 
IL S soles) Eek s(st& ) 230 5 (tw r(t.) 
7 -KA 1 KA 
a s(sto ) WE? : Sto yey a Ss(st oo) ES) 


and after taking the inverse Laplace transformation we arrive at 


x, faye : 3 f«r-1+e"*7) x, (kT)+(1-e7 %*)x, (kT)+ — rte" 7] 


x, ray2] SA ee” “Tx, (kt)te” Tx, (kt SE (dee) (Kr) 


A 
For x, (0) = x, (0) = 0.0 


myemeOr a Step input of ©(0)-= Y(T) “= OT )eseceece 4 ry 
we have 
= o2—7 eT 
Se - aT a -oxl iy 
ATGreee 7) A le (1+ «T)-1| 
aa. dal 
ho = ——— Ar hy = a NS: 
AT(C1=e ) MtiG@lae Ty 
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where xy C2) or: x, (2T) 20.0 


It is seen that a "delta'' term is present in the solution for the 
K's and the h's. This "delta" term arises from the fact that the 
error signal was attenuated by a 0.0121 factor before being sent to 


the A/D converter. This will be discussed in detail in Section 2.2.2. 


— 40 e => e a = e 
Now for the G(S) = areTay T = 0.3; where r = 5.0 we arrive at: 
K = 0.421 K, = -0.398 
O 1 
Hiee=.. 0.0252 h, = -0.0102 
O 1 


The h's computed were used in an earlier trial program. The computed 
K's were used in the final control program and as such were entered 
directly in the computer. The CDC-160 together with the CDC-168 were 
programmed to serve as the digital controller using these K. 8. The 
CDC-160 also operated as the sampler since the sampling period was also 
programmed into the computer. 

The real system response was also simulated using the same 
simulation scheme as Section 1.3. The simulation results are in 
pection 2.3. 

After controlling the real plant in a ''deadbeat response'' manner, 
it was decided to work on a time sharing routine for optimum control 
of two plants. The second plant was an analog simulation set up on 


the EAI TR-20 analog computer with a G(S) equal to The CDC- 


a 
s(stl) ~ 
160 was then used to provide optimum control to both systems. The 


design of the digital controller for the second plant followed the same 


theory used before. 


IL 





2.2 Development _of Sampled-Data Control Technique 


The step from the theoretical control solution to a working sampled- 
data controller involves the development of methods through which the 
requirements of the theoretical solution may be implemented in a real 
controller. This is a two-fold development involving the software of 
computer programming and the associated hardware requirements. Develop- 


ment in each area will be discussed in detail. 


2.2.1 Hardware Development 


Having developed, theoretically, the performance of the D(Z) 
portion of the control system block diagram, it was necessary to develop 
a hardware equivalent to the previously described mathematical model 
of D(Z). To accomplish this hardware development, the following equip- 
ment was used: 

(1) A two-channel digi eail to analog converter. (D/A) 
(2) A four-channel analog to digital converter. (A/D) 


(3) A Control Data Corporation 160 computer, in conjunction 
with a Control Data Corporation 168 arithmetic unit. 


(4) Operational amplifiers used for summing of signals coming 
in and out of the converter units. 


The equipment used was chosen, not because it was thought to be most 
efficient for the desired tasks, but because it was readily available 
in the Digital Control Laboratory of the U. S. Naval Postgraduate School. 
Comments on the relative merits of the equipment in performing these 
tasks will be included where appropriate. 

Development of the hardware equivalent of the mathematical model 


of D(Z) will now be discussed. It was decided immediately that the need 
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for analog hold circuits could be eliminated by replacing them with 
more accurate computer programmed holds. This is possible because of 
the operating characteristics of the converter units in conjunction 
with the 160 computer. First, the 160 computer can sample a converted 
digital input only when specifically commanded, and then it samples 
only one value. Similarly, the 160 computer can output a digital value 
for conversion only upon command, and furthermore, that value remains 
constant in the output register until it is changed by inserting a new 
value. Thus, by programming the computer to sample converted analog 
inputs and provide digital outputs for conversion at specified time 
intervals, hold circuits become unnecessary. In fact, the holds be- 
come an inherent part of the sampling rate of the computer and are 
varied automatically as the sampling rate changes. 

Examining Figure 6, one sees that the input to D(Z) is the error 
signal, e(t), which is the sum of the reference signal, r(t), and the 
negative control signal, -c(t). There are two methods by which e(t) 
may be transmitted to the computer. The first method is to sample r(t) 
on one channel of the A/D and to sample -c(t) on another channel. The 
two inputs could then be summed in the computer to obtain the error 
Signal. This method was discarded for two reasons. First, since it 
is impossible to sample r(t) and -c(t) simultaneously, the error signal 
computed would have some inherent error. Second, the A/D conversion 
time is of the order of 120 microseconds, and it was felt that this 
additional time delay might become prohibitive at sampling rates approach- 


ing 0.1 seconds. 


i 





The second method of supplying the error signal to the computer, 
and the one which was adopted, involved summing r(t) and -c(t) in an 
operational amplifier prior to transmission to the A/D. This method 
has the advantage of being fast, accurate, and readily adaptable to a 
restriction imposed by the physical characteristics of the converter 
unit. This restriction is that the converter units accept only negative 
voltages in the range of zero to minus ten volts. Thus, to handle both 
positive and negative values of e(t), the sum of r(t) and -c(t) is 
added by means of an operational amplifier to minus five volts. A 
zero error signal would then be sampled by the computer as an input of 
minus five volts, and an error signal of three volts would be sampled 
by the computer as an input of minus two volts. Similarly, the 
converted D/A value which is transmitted to the amplifiers of the 
system must be added to five volts prior to insertion in the system. 

A block diagram of the hardware system used to provide the 


equivalent to the mathematical model of D(Z) is shown below: 


| 160- 
> j= > = 


-c(t) 


To 


First 
Amplifier 


Figure 8 


h 


With the above hardware system established, the rest of the development 
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of a unit equivalent to the mathematical model of D(Z) must be 


accomplished by programming the computer. 


2.2.2 Software Development 


Programming the computer to accomplish the aforementioned task 

was done in a progressive fashion, starting with programs to accomplish 
the simplest control functions and expanding the basic programs to 
provide for accomplishing more complex tasks. This was done for two 
reasons. The primary reason was that at the outset the authors wanted 
to isolate quickly any faults in the digital control system, and to 
analyze carefully system response at each stage of control function 
complexity. A secondary reason for adopting this progressive program- 
ming technique was the fact that the 160 computer is rather limited in 
its sansa. Due to this limitation, the authors desired to bring 
the control capabilities of the equipment used to a maximum before 
computer saturation problems were encountered. Following is a discussion 
of these progressive programs. 

The first program written was intended primarily to check the 
operation of the hardware equivalent of the D(Z) under minimum capability 


requirements, and also to check the plant response with the D(Z) in the 


hohe 160 has limited storage capacity (4K), and a relatively slow 


memory cycle (6.4 microseconds). It can only multiply or divide in 
conjunction with the 168 and the 900 microsecond time for such operations 
is quite slow. 


Jie 





loop. The essence of Program 1 is to enable the plant to run as a 
continuous system with the D(Z) block in the loop. Mathematically, 
such a system may be described as D(Z) = 1., and the sampling period, 
T, approaching zero (approximately 200 microseconds). 

The programming aspects of Program 1 were quite simple. The 
program was an iterative loop in which the current error signal was 
sampled from the A/D converter, read into the computer, stored 
temporarily, and output to the D/A converter for return to the system. 
A flow chart for Program 1 is shown below. A complete text of Program 


l is presented in Appendix VI. 


Program 1 


from Temp 
Storage 





Figure 9 
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The results of the test using Program 1 were quite satisfactory. 
System responses were obtained from step, ramp, and sinusoidal inputs, 
and these responses matched those obtained from the system when operated 
as a purely continuous system with unity feedback. This, of course, was 
predictable with D(Z) = 1., and T approaching zero. However, the results 
do validate the selection of hardware used to implement the realization 
of the mathematical D(Z). 

The results of Program 1 opened another avenue of approach to the 
problem of implementing the variable gain samplifier method of digital 
control. It was noted, when working with the system in the Digital 
Control Laboratory, that the system was subject to many random inputs 
of small magnitude from various sources of noise. Program 1, when 
tested, maintained stability, and accurate and fast response in the 
presence of these noise sources. It was further noted that the advantages 
of the variable gain amplifier method in regard to minimum time response 
are lessened when the inputs are small and random. 

In the light of the foregoing facts, a temporary control philosophy 
was developed. This philosophy states that for step inputs which exceed 
in magnitude a threshold, a standard variable gain amplifier control 
solution would be employed. For small inputs, such as noise, the system 
would be controlled as it was with Program 1. The threshold would be 
set at a level which would prevent the implementation of a solution of 
the type described in Section 1.2 for the random noise inputs, which 


were restricted to a relatively low magnitude. In other words, the 
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system would run in an essentially continuous mode D(Z)=1., until 
receipt of a valid step input, at which time it would be controlled 
by a variable gain amplifier digital solution. 

To implement this type of control philosophy, an intermediate 
program, Program 2, was written to modify Program 1 to provide for 
detection of a definite step input and an exit to a solution. This 
detection is accomplished by comparing the new value of the sampled 
error signal with the immediately preceding value of the sampled error 
signal. If the difference between these values is greater than the 
established threshold, a step input is detected. To assist the reader 
in understanding the programming procedure for determining the sign of 
a detected step input, the following table of octal number values for 


analog voltages into the converter is presented: 


ANALOG OCTAL 

VOLTAGE NUMBER 
0 4000 
-1 4632 
-2 5463 
=9 6314 
~4 7144 
-5 0000 
-6 0631 
ai 1463 
-8 2314 
-9 3144 
-10 3777 


A flow chart for Program 2 is presented on the following page. A 
complete text of Program 2 is presented in Appendix VI. 
The results of Program 2 were excellent. It proved to be a fast 


and accurate method of sensing step inputs having a magnitude greater 
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than the threshold level. In one sense Program 2 was too sensitive. 

When certain knife switches were used to provide zero to five volt step 
inputs, negative transients appeared at the instant the switch was thrown. 
Program 2 was so sensitive that it sensed the transient as a negative 
step input. This imposed the requirement that step inputs be obtained 
from an electronic switching device or from a high quality knife switch 
which eliminates the aforementioned transient. 

The next step in the development of the temporary control philosophy 
involved utilizing the computer in a real time mode to implement a 
variable gain amplifier solution. The iterative loop of Program 2, 
which maintains continuous control until a step input is sensed, employs 
the computer in a free-running mode. The variable gain amplifier solution, 
however, requires that the computer insert control voltages of definite 
real time length into the system. To do this, the computer must output 
a control voltage to the system, hold the voltage for a specified length 
of time, and then output another control voltage. Program 3 was written 
to develop techniques for accomplishing this task. 

Program 3 was a simple program which called for the computer to 
output a given voltage, delay through a timing chain, and output another 
voltage. Theoretically, the length of the delay could be calculated 
from the execution time for the instructions which constitute the timing 
chain. However, exact timing data for D/A conversion were not available. 
Therefore, approximate timing chain delay was calculated and adjusted 
by experimental results for exactness. The timing chain was formed by 
constructing an iterative loop in which an index was increased by one 


on each pass and compared to a preset total. When the index equaled the 
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preset total, the time delay was complete. A flow chart for Program 3 
is shown on the following page. A complete text of Program 3 is 
presented in Appendix VI. 

In testing Program 3, it was decided to design a basic delay block 
of 0.1 second. With an accurate delay block established, delays which 
were integer multiples of 0.1 second could be achieved by iterating 
through the basic delay the required number of times. Delays of less 
than 0.1 second could be achieved by setting the preset total at a proper 
fraction of that required for a 0.1 second delay. 

Satisfactory results were obtained from Program 3 in that a delay 
of exactly 0.1 second was produced between the output of the first and 
second voltages. For this delay, the preset total was set at 3403... 
Results were checked on a high speed Mark II Brush Recorder and an 
oscilloscope. 

Having developed Programs 2 and 3, the basic tools for implementing 
a solution within the constraints of the temporary control philosophy 
were ready. By using Program 2 to maintain continuous type control in 
the absence of an input signal and to sense a step input, and by employing 
the features of Program 3 to output the required controlling voltages, 
the desired solution to a step input may be attained. To test the basic 
validity of their approach, the authors wrote Program 4 as a first trial 
of the variable gain amplifier control method. 

Program 4 was an unsophisticated approach which started with the 
system being controlled in the undisturbed state by Program 2 which will 


hereafter be referred to as the Sense Loop. When an input is sensed, ho 
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is transmitted immediately to the system. ho is held for one sampling 


period, and then h, is transmitted and held for one sampling period. At 


1 


the completion of the second sampling period, the system returns to 





continuous operation as the computer returns to the Sense Loop. To 

insure a return to the continuous mode of operation at the completion 

of the second sampling period, it was necessary to ''bootstrap" the Sense 
Loop by sensing the current error and inserting it in the Sense Loop as 
the "old value". Were this procedure not followed, an imperfect solution, 
one which did not have zero error at the end of two periods, would cause 
the Sense Loop to sense another step input and thus cause the system 

to go into unstable, erratic operation. 

The unsophisticated aspect of Program 4 stems from the fact that 
the h's were not computed by multiplying e(t) by K during the solution, 
but were pre-calculated and inserted in memory as constants to be trans- 
mitted to the system at the proper time. K, and K, have been determined 


0 1 


in Section 2.1. ho was simply Ky times en (t), and h, was K, times e, (t) 
(predicted). Due to the fact that r(t) and -c(t) were transmitted to 
the summing operational amplifier after passing through dropping resistors 
in the front end of the system, the magnitude of the error signal was 
considerably decreased. For a five volt step input, ep (t) sensed by 
the computer was 0.0600 volts and predicted e, (t) was 0.0258 volts. On 
the basis of these values, ho = 0.0252 volts and hy = -0.0102 volts. A 
O.3 second sampling interval was used. 

A flow chart for Program 4 is shown on the following page. A 
complete text of Program 4 is presented in Appendix VI. 

Complete results of the responses obtained when using Program 4 


will not be presented because these tests were in the nature of feasibil- 


ity checks. The results may be summarized by stating that they provided 
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definite proof of the validity of the variable gain amplifier approach, 
and led the authors to adopt certain techniques which were used in 
later programs to improve performance. Some of these lessons learned 
from these tests will now be discussed. 

Accuracy problems were caused by the small magnitude of the error 
signal received at the D(Z) portion of the system for a five volt 
step input. Due to this fact, the magnitude of h's to be transmitted 
to the first amplifier was so small that the resolution of the D/A 
converter prevented sufficient accuracy for these h's. The resolution 
for the converter is 0.0024 volts per octal number. With h's of the 
magnitude of 0.02 volts, accuracy was limited to the first digit of 
the h. This condition was barely satisfactory due to the stringent 
demands of the variable gain amplifier method for accuracy in the h's. 
This demand was demonstrated when solutions were simulated on a 
general purpose digital computer. This demand was discussed more 
Broroughly in Section 1.3. 

To correct the problems associated with converter resolution, 

the obvious solution was to increase the voltage levels in and out 
of the converter. The simplest method of accomplishing this would 
have been to eliminate the dropping resistors before the summing 
amplifiers. However, when this was tried, system noiSe increased to 
a level which was not tolerable. The method of achieving increased 
accuracy, which produced the best results, was to multiply the error 


signal by ten in the summing operational amplifiers prior to A/D 


Fis, 





conversion, and to divide the D/A converted voltage by ten in the 
summing operational amplifiers before transmittal to the system. This 
enabled the use of h's of ten times computer value for conversion 
which alleviated the resolution problem. This scheme also maintains 
operation in the continuous mode of the Sense Loop to continue as 
before. For all further tests, the hardware shown in Figure 8 was 
modified to provide for the proper multiplication and division in 

the summing operational amplifiers. 

A second lesson learned in this test program related to the 
"bootstrap'' procedure discussed previously. This feature was 
originally inserted in the program to prevent system runaway for 
incorrect solutions resulting from tests using improper h's. It was 
found, however, that this "bootstrapping" return to continuous operation 
at the completion of the second period was an integral part of a 
satisfactory solution. This is due to the fact that a real system 
has inherent non-linearities such as coulomb friction and backlash. 

In view of these non-linearities, a perfect theoretical solution can 
not be obtained. (It is possible to approach more closely the perfect 
solution by slight modification of the size of the h's.) Thus, a 
rapid, smooth return to the continuous mode after the second sampling 
period is essential to a satisfactory solution. It is noted that 

by properly adjusting the size of the h's, the error at the end of 

the second sampling period is small (less than 10% of the step size), 
and the use of the continuous mode reduces error to zero rapidly. 


Specific examples of the foregoing may be noted in the presentation 
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of results in Section 2.3. At that time, the reader may note that 
the necessity to go to the continuous mode to bring error to zero 
does not significantly degrade system performance in comparison to a 
theoretical solution. 

With the information derived from Programs 1 through 4, the 
authors developed Program 5 which was intended to culminate efforts 
in the development of the temporary control philosophy. Program 5 was 
a generalized version of Program 4, the primary difference being that 
in the newest program the h's were computed during the solution 
instead of using precomputed values. The solution technique may be 
summarized by stating that each h is computed by sampling current 
error and multiplying by the appropriate gain constant (K,) by using 
the CDC-168 arithmetic unit in the multiply integer mode in conjunction 
with the CDC-160 computer. 

Programming arithmetic operations for these two units is quite 
straightforward when Subroutine Arith, written by Professor M, L, 
Cotton, is used. However, the characteristics of the CDC-168 required 
manipulation of the sensed error signals. Specifically, for an 
arithmetic operation, the 168 requires two 22-bit operands, and 
supplies a 22-bit solution. Each operand is composed of two 12-bit 
160 words, with the least significant half of the number in an even 
numbered cell and the most significant half of the number in its odd 
numbered mate. The first bit in each cell pair is a sign bit. It is 
this sign bit in the first bit of each of the two words composing the 
operands that require the manipulation. For the K multipliers which 


remain constant, proper values may be inserted into both storage cells 


ol 





manually before the program runs. However, for the sensed error 
signal, negative numbers received from the converter create a problem. 
For the 168 to function properly, the sign bits for both portions of 
the operands must agree. To handle this problem, the following proce- 
dure was employed. To derive the proper h, the product of K. and e(t) 
is required. After e(t) is sensed, it is placed in cell 0024 as the 
least significant half of an operand. The most significant half, 

cell 0025, is preset to zero. For negative step inputs, the error 
signals read in from the A/D converter are positive numbers, and the 
multiplication process may proceed immediately. For positive step 
inputs, however, the error signals sensed from the converter are all 
negative numbers. To maintain sign bit consistency with cell 0025, 
this negative number is complemented and placed in cell 0024. The 
multiplication operation then proceeds. The product, of course, is 

of the wrong sign so it must be complemented prior to being transmitted 
to the system. This procedure will be defined completely in the flow 
chart for Program 5. 

The other major innovations of Program 5 involved attempts to 
obtain greater accuracy in the solution and maintain high resolution 
in the converter. As previously discussed, the multiplication and 
division by ten in the operational amplifiers were included. To 
achieve greater accuracy in computing the h's, 100 times K, was set 
in the computer as the multiplying constant. The products were then 


divided by 100 prior to being output to the D/A. 
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As before, a sampling interval of T = 0.3 seconds was used. 


The gain constants for this system were computed to be: K,=0.421 and 


0 
Ky = -0.398. A flow chart for Program 5 is shown below. A complete 
text of Program 5 is presented in Appendix VI. 
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As noted before, Program 5 marked the culmination of efforts 
in conjunction with the temporary control philosophy. Although the 
results of Program 5 were good, and gave valuable insight into digital 
control technique, Program 5 was merely a vehicle through which the 
validity of variable gain amplifier theory could be tested. The 
temporary control philosophy, although perfectly satisfactory for 
control of a single system, is unsatisfactory as a general philosophy 
for digital control because it excludes utilization of one of the 
primary advantages claimed for digital controllers: simultaneous 
control of many systems by the same computer. The authors believe 
that the realistic application of digital control methods is closely 
tied to the capability for simultaneous control of many systems. 
Toward this end, a final control philosophy was developed. 

The goal of this final control philosophy was to achieve this 
simultaneous control. This goal implies that the computer be "time- 
shared" among the systems to be controlled. Therefore, the significant 
difference between the temporary and final control philosophies is 
one which enables the computer to time-share its control function. 

Before discussing the final control philosophy in detail, 
specific performance goals for this philosophy will be discussed. 

The authors wanted to develop techniques for controlling two or more 
systems simultaneously by employing one digital computer on a time- 
share basis. The control was to be achieved by the variable gain 


amplifier method. To achieve generality, it was desired to include 
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the capability for the systems to be controlled at different 

sampling rates. Furthermore, it was considered necessary to maintain 
the capability to handle simultaneous step inputs to each system 
without performance degradation. 

It may be recalled that the temporary control philosophy called 
for a return to a free-running "continuous" mode at the completion of 
the solution. The major difference in the final control philosophy 
is that at the completion of the second sampling period, the computer 
returns to a mode in which D(Z) is still 1.0, but instead of a sampling 
interval approaching zero, a Sampling period of definite finite 
length is employed. Although this innovation does not change the 
response characteristics of the control solution, it marks the change 
from a partially discrete, partially continuous system, to one which 
is completely discrete. It is this change which enables the develop- 
ment of the required time-sharing computer operation for multiple 
control. 

Having discussed the major difference in the control philosophies, 
the final control philosophy will be covered in detail. As before, 
with no signal input to either system, the computer resides in a 
sense Loop. In the Sense Loop each system under the control of the 
computer is sampled and tested for a step input. If a step input 
greater than a threshold is not detected, the exact signal sampled 
is returned to the system. (D(Z) = 1.0) Each of the systems under 
control is handled in this fashion. When all systems have been 


sampled, the computer goes through a basic delay. This basic delay 
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is so named because it is a building block for all delays required 
for implementing variable gain amplifier solutions of any length on 
any system controlled by the computer. All required solution delays 
are formed as integer multiples of the basic delay. It is noted that 
the basic delay is much longer than the time required to sample and 
test all systems being controlled. 

When a step input is received by one of the systems, an error 
flag for the appropriate system is set and the computer exits. from 
the Sense Loop to a Solution Routine. In the Solution Routine, the 
first control voltage for the flagged system is transmitted from the 
computer, and all other system flags are checked for possible simul- 
taneous step inputs. If not, the computer goes through one basic 
delay and then senses all other systems for possible step inputs. If 
no other systems have received signals, the computer returns to the 
start of the Solution Routine and begins to count the number of times 
the basic delay has been entered. After each pass through the basic 
delay, all other systems are checked for possible step inputs. When 
the pass count reaches a number which indicates that the first control 
voltage has been in the proper length of time, the computer computes 
and outputs the second control voltage, and a new pass count begins. 
If, at any time during the solution, another system ‘receives a step 
input,it is flagged, and on the next pass through the Solution Routine, 
a solution similar to the one explained above is begun. 

When a variable gain amplifier solution for one of the systems 


- 


is completed, the computer jumps to a Return Routine. In this Return 
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Routine all counters for the appropriate system are zeroed and the 
error flag is set to zero. The system is then bootstrapped by the 
procedure previously explained, and returned either to the Sense Loop 
or the Solution Routine, if the error flag of another system is set. 
In any case, after being bootstrapped, the system returns to a mode 
in which the characteristics of the control function are D(Z) = 1.0, 
and T = the basic delay. 

The characteristics of the final control philosophy may be 
summarized by stating that any single system which is being controlled 
operates in a state in which D(Z) = 1.0, and T = the basic delay when 
there is no step input above the threshold level. Thus, for small 
step inputs and noise inputs, the system would be controlled in the 
above mode. When a step input greater than the threshold level is 
received, a variable gain amplifier solution is executed and the 
control mode then returns to that specified above. In the time 
sharing mode, each system controlled is operated in this same manner, 
with the option of making the solution sampling period any integer 
multiple of the basic delay. 

To test the final control philosophy, Program 6 was written. 

For purposes of this test, two systems were time-shared. One of the 
systems used was the real system used in previous tests. The second 


system was an analog simulation of a plant having the following 


1 


Try The basic delay was established 


characteristics: G(S) = 
as 0.025 seconds. The solution sampling rate for the real system 


was 0.3 seconds; for the analog system, 1.0 seconds. The only 


Sh) 








restriction for Program 6 was that it was written to control only 

two systems simultaneously. The basic delay time may be chosen 
arbitrarily, and the solution sampling rates may be any integer 
multiple of the basic delay. Computation of ho and hy was accomplished 
in the manner of Program 5. The hardware arrangement for the real 
system was the same as that employed in conjunction with Program 5, 

and the summing amplifier arrangement in and out of the converter for 
the analog system is identical to that of the real system. 

A flow chart for Program 6 is shown on the following pages. It 
is broken down into three sections: Sense Loop, Solution Routine, 
and Return Routine. A complete text of Program 6 is presented in 
Appendix VI. 

To simplify the flow chart for the Solution Routine, the negative 
step solution for System 1, and the positive step solution for System 
2 have been omitted. The omitted processes are identical in form to 
those shown. Furthermore, the "D/A hy," blocks, shown in the flow 
chart as predefined processes are handled exactly as in Program 5. 

Before presenting the flow chart for the Solution Routine, the 


following indices are defined: 


I = the number of passes through the basic delay for the 
current h for System l. 


J = the number of h's transmitted to System 1. 
P = the required number of passes through the basic delay to 
achieve the specified variable gain amplifier solution 


sampling rate for System l. 


K = the number of passes through the basic delay for the 
current h for System 2. 
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iL. = the number of h Ss tEtansmijtmedero mySeem 2. 


Q = the required number of passes through the basic delay to 
achieve the specified variable gain amplifier solution 
sampling rate for System 2. 


The flow chart for the Solution Routine is shown on the following 


page. 
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The Return Routine is shown on the following page. The Return 
Routine will be shown for System 1 only. That for System 2 is identical 


except for the fact that different indices are zeroed. 
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2.3 Results 

The results of Program 6 were excellent. Both systems were controlled 
with predicted speed and accuracy. The real system response was equivalent 
to that obtained from utilization of Program 5. It is interesting to 
note that the system simulated on the analog computer exhibited exact 
theoretical response to the variable gain amplifier solution. This, of 
course, is due to the fact that there were no non-linearities in the 
Simulated system. However, the exact theoretical solution obtained does 
validate the final control philosophy employed. 

The time-share program was tested for all possible time combinations 
of step inputs to the two systems. In no case was system response de- 
graded once the solution was started. In the worst case, there was a 
25 millisecond delay between a step input and the start of a solution. 

This occurred when a step was entered just as the computer started a 

pass through the basic delay in the Sense Loop. The delay, in this worst 
case, which is a random occurrence, was not discernable on the Mark II 
Brush Recorder used to measure response. 

On the following pages, the response curves for the two systems 
controlled in the time-share mode by Program 6 are shown. Also shown 


are the theoretical response curves for the real system. 


2.4 Extension of Results 

Due to the limitations imposed by the D/A, only two systems could 
be simultaneously controlled in the Digital Control Laboratory. To 
fully substantiate the validity of the final control philosophy, one 


must consider the feasibility of this philosophy in controlling more 
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Figure 20 
Real System 


Upper Graph Velocity vs. Time; Lower Graph Position vs. Time 
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‘Figure ale 
Analog Simulation 


Upper Graph Velocity vs. Time; Lower Graph Position vs. Time 
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than two systems. Let us consider this feasibility with respect to the 
equipment used throughout the laboratory tests, assuming only that the 
number of channels in the converter may be arbitrarily expanded. (The 
authors feel that the problems associated with more difficult tasks in 
the computer field are too often dismissed with a bland inference that a 
bigger, faster computer can do any job. The dollar economics of the 
control problem preclude this approach.) 

Using Program 6 as a basis for time-sharing control of many systems, 
the CDC-160 has sufficient memory capability to handle 12 systems 
simultaneously. By sub-routinizing wherever possible, this number could 
probably be increased to 15. For 15 systems, however, the sampling 
time, governed primarily by the analog to digital conversion time of 
120 microseconds, becomes an appreciable proportion of the basic delay 
time. In the Sense Loop, for example, the time to sample A/D, check for 
step inputs, and return D/A would be approximately two milliseconds 
compared to the 25 milliseconds of the basic delay. This added delay 
in the Solution Routine would produce significant errors in the variable 
gain amplifier solution. If the basic delay could be satisfactorily 
increased to 50 or 100 milliseconds, the two milliseconds required to 
sample the 15 systems would become insignificant and not seriously 
degrade the solution. The authors feel that employment of a 50 or 100 
millisecond basic delay, which limits the fastest solution to a step 
input to 100 or 200 milliseconds, is not unreasonable in most practical 
applications. 

The foregoing discussion has made one tacit assumption. That 


assumption is that the number of simultaneous inputs is limited to three 
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of the 15 systems. This is due to the 900 microsecond multiply time 
required by the CDC-168. This, of course, excludes the desired 
generality for handling any step input at any time for any system. 
With the equipment used, the 160 in conjunction with the 168, there 
is no way to compensate for this problem. 

At this point it is worthwhile to look at the advantages of using a 
different computer. Specifically, the CDC-160A would adequately solve 
the problem mentioned above. The 160A is a computer quite similar to 
the 160 in all respects with the additional capability for fast 
multiplication and division. Economically, the 160A is comparable in 
cost to a 160-168 tandem, and makes a much more compact unit. 

On the basis of the foregoing tests and theoretical extension, 
the authors feel that with a small, relatively inexpensive computer, 
such as the CDC-160A, it is quite feasible to control up to 15 systems 
Simultaneously. The control capability will vary with the capabilities 
of the brand and type computer selected. The main thesis, however, 
is that an excellent, multi-system, digital control scheme may be 


implemented without going to the larger general purpose computer. 
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CHAPTER IIl 


3.L Conclusion 


It was noted in the Introduction that the basic goals of this 
thesis were the controlling of hardware by digital methods, and the 
developing of techniques for implementing time-sharing control of many 
systems from the same computer. The results obtained have been detailed. 
In this section the most important findings will be summarized and their 
relevance to the general feasibility of digital control systems will 
be noted. 

The most significant aspect of controlling the hardware was the 
necessity for going to the "continuous" mode (D(Z) = 1. and a small T) 
to bring the error to zero at the conclusion of the solution to a step 
input. This characteristic, caused by the non-linearities of the system, 
would seem to be a significant disadvantage to this type control. 
However, this was not the case. In the worst case, the servo position 
was within ten percent of the desired final position at the moment of 
the switch to the "continuous" mode. Furthermore, at the switch time, 
the servo velocity was in the proper direction at a decreasing magnitude. 
This condition permitted the system to settle rapidly to the desired 
position. 

Of course, to realistically evaluate this control method, one 
must establish a performance criterion as the basis for comparison with 
other control methods. The criterion selected was one which measured 


minimum time to reach a position within ten percent of final value and 
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to remain within this ten percent boundary while settling. On the 
basis of this criterion, the time for a unity feedback continuous 
system to meet the conditions was 1.2 seconds, and 800 milliseconds 
with the addition of optimum tachometer feedback. These figures were 
obtained with the amplifier gain potentiometers set at the same levels 
used for the digital operation. By increasing the amplifier gains, a 
time of 500 milliseconds was obtained with tachometer feedback. Using 
a variable gain amplifier solution sampling rate of 300 milliseconds, 
a time of 600 milliseconds was required to meet the criterion. The 
relatively low saturation level of the amplifiers precluded the use of 
sampling rates of 100 or 200 milliseconds, but on the basis of responses 
obtained at rates of 300, 400, and 500 milliseconds, it follows that 
improved amplifiers would allow a solution time of 200 milliseconds. 
In general, the fastest speed of response obtainable by the variable 
gain amplifier method is limited by the saturation level of the 
amplifiers. Similarly, the response speed obtainable by a bang-bang 
technique is limited by this saturation level. 

On the basis of these results, the performance characteristics 
of the digital control method indicate that this method is competitive 
with other control methods. With improved amplifiers, response speeds 
from the digital methods can be much faster than those obtainable from 
continuous methods. Optimum bang-bang techniques approach the response 
speeds obtainable with the variable gain amplifier method. 


There are, of course, many criteria by which control systems 


5). 


may be judged. The criterion above is one of speed of response. 
Although other criteria were not evaluated, it is noted that the 
digital technique is, in general, quite adaptable. For example, the 
basic hardware and software employed by the authors is readily adapt- 
able to a criterion which calls for minimum fuel expenditure. 

The results of the time-sharing program have been fully discussed 
and the feasibility of extending the program to control many systems 
has been outlined. The basic approach of utilizing a small computer 
has precluded a discussion of a further extension of the technique. 
That is, using a large, general purpose computer as the digital 
controller, with other computing tasks being time-shared with the 
control function. This is feasible with no change in the basic program 
heretofore used in conjunction with the CDC-160. In examining this 
program, one can readily see that the vast majority of time is spent 
in the basic delay block. Time spent in the basic delay is ideal for 
carrying out other computing tasks. This could readily be implemented 
on a general purpose computer with a real time clock and interrupt 
capability, and an executive routine to direct the proper sequence of 
computing tasks not associated with the control task. The ramifications 
of this large computer capability are indicated by Slaughter and 
Lackowski in a paper presented to the 1963 National Convention on 
Military Electronics in which they state: "In those cases in which a 
digital computer is available, digital control requires less hardware 


than does conventional control."| 4 | The presumption accompanying this 
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statement was that the digital computer was required for other tasks, 
and the controlling tasks could be satisfactorily time shared. 

In view of the performance characteristics of the variable gain 
amplifier method and the time-sharing capabilities developed, the 
authors have come to the conclusion that for many applications, digital 
control methods are presently feasible and superior to other control 
methods. Due to the present cost of computers, these applications are 
presently limited to large processes or areas where the advantages of 
digital control are overwhelming. A specific example of such an 
application is in large processes such as those in the chemical or 
petroleum industries in which multiple systems must be simultaneously 
controlled. Similarly, a single large process with multiple inputs 
may be controlled by digital methods. A military application would 
be the extension of the NTDS to include the servo control of gun 
mounts and missile launchers from the master computer. A final 
application might be in the aerospace field in which the new techniques 
for optimum control on the basis of minimizing a given cost function 
are quite adaptable to digital techniques. 

Although immediate feasibility for digital control is limited 
to processes of the type mentioned above, the decreasing size and 
cost of digital computers increase their field of application. 

The investigation for this thesis was limited in scope to permit 
the extensive study of a single method. The results from this study 
may, in general, be extended to other digital control techniques. The 
authors feel that the results of this study of digital control validate 
its basic feasibility, and indicate the value of further study in the 


field. 
5/ 
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APPENDIX I 


DESIGN OF DIGITAL CONTROLLERS AND SIMULATION RESULTS 


In addition to the design of a digital controller for a step input 


Ae 
s(s+1) 


designed. Referring to Figure 1 and Figure 2 in Section 1.2 we arrive 


to a plant, a digital controller for a unit ramp input was also 


at the following equations. 
T = 1.0 


x [cee | = 1-0. 368K, | x, (kT) + 0.632x, (kT) + (0.368K,)x(kT) 
By [+17 | = (-0.632K,)x, (kT) + 0.368x, (kT) + (0.632K, )r(kT) 


For a unit ramp input r(kT) = k. 


k = 0; x, (0) = x» (0) = 0.0 
x, (T) = 0; x, (T) = 0 
k= 1 
x, (2T) = 0.368K, 
x, (2T) = 0.632K, 
To obtain a solution we must set x, (nT)=x (nT) and x, (nT)=1. We there- 


fore must proceed to the next period in order to solve for the given 
conditions. 


Therefore k = 2 


x, (3T)= [ 1-0.368K, (0. 368K, ) ae (0.632)°K, ar 2 (0. 368K, ) 2 3.0 


Z 


x, (3T)= | -0.632K, | (0.368K,) + (0.368) (0.632K, )+2(0.632K,) 41.0 
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from which 


Ky = 3.328 K, = 0 peul 


and 


h(T) = 3.82 ; h(2T) = 0.183 


It should be noted that for this system to maintain a zero error after 


it has arrived "home'' a constant output of h(nT)=1 must be fed into the 


continuous system. This will keep the system moving at a unit velocity 


and thus maintain a zero system error. 


A digital controller using the same design theory was developed 


1 


woe (s+) (st2) 


The derivations and explanations follow. 


Same as that in Section 1.2. 


i 


Plant: G(S) = CD CE 


Sampling Period: T = 1.0 


Block Diagram: 





(or wero 
; | | 
a Gs (T) iL | D(S) 
+ Ak | T | 
| | 
Nes weet Fone eerie 
D(Z) 
Figure 24 


State Transition Flow Graph: 


_ 


: . K(t,) z S 5 
re) e(t”) h(t?) H(s) 


Figure 25 
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1 


(s+1) (s+2) 


The procedure used is the 


plant in response to both unit step and ramp inputs. 


x, =C(t) 





xy | ct) ] 
Xo | ct1)2| 


2, te Ta (5-077. 50777)| x, (KT)+(e7*~e7°7) x, (KE)+ 
K,(.5-e + se Ce) 


j2e7 7-267 Ti Ce" 
K,(e7-e **) r (kT) 


Al 


eT) | x, (et) # (207-0) x, (kt )+ 


For T = 1.0 


ey [oone] 


Xo [ erie | = |-.466 - 233K, x, (kT) = -098x, (kT) or - 233K, r (kT) 


|.601-.1995K, | x, (kT) ai -233x, (kT) I -19995K,. vel AB 


(L) Let r(kT) = 1.0 


k=0 ; x, (0) = x, (0) = 0 
x, (T) = -1995K. 5 Xo (T) = -233K. 


These two equations clearly cannot be solved for the final conditions 


Of: x, (T) = a0: x, (T) = 0.0 


Therefore let k = 1 


Ite 


x, (2T) 0.1742K - 0.0398K K + 0,1995K 


1 l ae 


I> 


x,(2T) = -0.11575K - 0.0465K K, + 0.233K 


1 1 


Prem which: 
fer = 5.65 ke = 6.73 
O 1 


and 


h(O) = 3.65 ; h(T) = 1.832 
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Note: 


steady-state error to a step input. 


Since the plant involved is of type zero it will have an inherent 


To eliminate this error a "steady 


state h'' will have to be put into the system after the error has been 


initially reduced to zero. 


value theorem was used as follows: 





Assume: h(t) = at + b 
a b 
HG = a + me 
(oie — a 
=| (st1) (s+2) 
d_ | (bsta)e>" 
a x, (oo ) ds ? 
s +3s+2 
2b+2at-3a 4a 
CO ee = 1.0 
or a=0; b=2 
Poe (ih) n = 2,3 
(II) Now let r(kT)=k 1..e. unit ramp 


In this case the final values desired are: 


for n the smallest possible integer. 
FOr: ik.=.0; x, (0) =x, (0) =0 

x, (T)=0 2 x, (T) = 0 
Therefore let k=1 


x, (2T)= 0.199K, 
x, (2T)= 0.235, 
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To determine this "steady state h" the final 


x, (nT) =k; x, (nT)=1.0 





These two equations cannot be solved so we must proceed to k=2. 


= 
x, GT)= 0.17425K, - 0.0398K, K, ate 0.399K, = 3.0 


4 
x, 3T)= -O.11575K, - 0.0465K)K, ate 0.466K, = 1.0 


From which 


Kaw = 7.85.2 Kare F671 
and 


RCL) .= 7.8575 hry = Gale 


Again to keep the output equal to the input for this type zero system 
we must output a "steady state h''. From the final value theorem this 
AS: 

hQnd) y=) 2e-oo on Bee. ok: 


The following pages contain the graphical simulation results for 


a 
s(s+1) 


responses of a type zero system (G(S) = 


the type one system (G(S) = ) to a unit ramp input and the 


1 
(aD ) to both unit step 


and unit ramp inputs. 
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vane 
CEE 


—RAOSCALE + 4.0QE+A2 UNITS“ INCH, 
M-SCALE = 1.28£+89 UNITS/INCH, 


BROWNE DIGIT 12 © <aey t=} 
RUN 4 OUTPUT US. TIME 


Figure 26 
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S 

a 
MR 








A-SCALR = J,ARE+R80 UNITS“ INCH 
M-SCALIC = 0,80 -81 UN] TS7 INCH. 


RBKOWNE DIGIT IR 
KUN J 


VELOCITY US.TIME 


Figure 2/7 
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A-SCAL = 2,QQRE-BJ UNI TS~ INCH, 
“TSCA = ORR -21 UNITS“ INCH, 


BKOWNE DIGIT 1K 
KUN J EDOT US. ERROR 


Figure 28 
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G10 
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BQN. QQ5 


A-SCALE = 1,Q0£+20 UNITS/INCH, 


M~SCALE = 2.Q8E-24 UNI) S~ INCH, 


BROWNE DIGIT 2 
KUN J 


1 
G(S) =(serytsrzy or (kT)= 1.0 


OUTPUT Vs. TIME 


Figure 29 
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QQ3 BOY Q25 


K-SCALE ~ 1.Q8E+82 UNITS“INCH. 
4-SCALE = 2.Q86-241 UNITS” INCH, 


BROWNE DIGIT 2 
KUN 4 


VELOCITY US. TIME 


Figure 30 
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A-SCALE = 2.828E-Q4 UNITS INCH, 
M-SCALE = 2.8Q6°21 UNITS/ INCH, 


BROWNE DIOCIT 2 ; 
KUN J EDOT US. ERROR 


Figure 31 
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X-SCALIL = 1,.Q2E+R2 UNITS” INCH. 
M-SCAL = 4.208422 UNITS/INCH. 


Seon OICIP CR “Ose een “=. 
KLIN J OUTPUT US. TIME 


Figure 32 
FAO, 








A-SCALR = 1,Q2£4+28 UNITS/“INCH 
M-SCALK = SRABE-21 UNITS“ INCM 


BROWNE DIGIT 2R : 
KUN J VELOCITY US.TIME 


Figure 33 
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| ~ / 
| PKL 


A-SCALI = 2,026 -2) UNITSINCH. 
“4-SCAL = S2QQE-B1 UNITSZINCH. 


BROWNE DIGIT 2R | 
KUN J CDOT US. ERROR 


Figure 34 
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APPENDIX IT 
SIMULATION PROGRAM 
The following pages contain the simulation program used in 
verifying the digital controller design for a "deadbeat response" to a 


unit step input. The programs for simulating a unit ramp response are 


Similar. 
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MeL se BROWNE Esko. 


10 
tal 
eZ 


i> 
14 


PROGRAMs © 10G sale 
DIMENSION X(30) sxXDOT (40 eal >) 
Ct 103 = 15.0 

CALL INTEGI S(T sX>sX DOS) 
INPUT=1e0 

LF A=) eee 1) )) LO ieee 
HOLD® C2) 

GO TO 14 

HOR = 986 (3) 

GO TO 14 

lIF (7=20*GG 1) 1) sie 
HOLD = C(4) 
XDOT(2)=HOLD-X( 2) 
XDOT(1)=X(2) 
ERROR=1¢0-X(1) 


EDOT=-X(2) 
X( 3) =ERKOR 
X(4)=EDOT 
C(11)=50.0 
GO. [Ca 
END 

END 
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BROWNE DIGIT 1 
ONE RUN IS CALLED FOR | 


INPUT DATA RECORD 


GROER OF EQUATIONS = 2 
INITIAL TIME = eOQO00E+C0 
FINAL TIME = .5000E+CO 
STEP SIZE = .2000E-C3 
THE NCN-ZERO CONSTANTS, CCI), ARE 
C( 31) = .1000E+CO 
Cl. 2) =i OS secs 
C( 3) = -.9580E+C2 
ALL THE INITIAL CONDITIONS ARE ZERO 
THE COLUMN HEADINGS AND THE CCRRESPONDING VARIABLES ARE 
TIME X{ 0) 
OUTPUT X( 7) 
VELOCITY X{ 2) 
ERROR X( 3) 


THE INDIVICUAL GRAPH TITLES AND THE CORRESPONDING VARIABLES ARE 


OUTPUT VS. TIME X( 1) VS. Xf Q) 
VECECTE evs. [lhe X( 2) VS. Xt CO) 
ECOT VS. ERROR XC 4) VS. XC 3) 


iD 





BROWNE DIGIT 1 


VELOCITY ERROR 


OUT PUT 


TeEMe 


m™DODOVI0D00000 
ODOdOVIO00000 
+-eeee ee + + + 
LULU LU bay WL Lo Lt 
Ost Nit here oop 
OP- ODP Oe OAD 
OtRAOOMEO st 
OOP Ute Pe rm th 
mm OOOO OP Ou 


ee ,e0848¢e 6 @ 


OC en en mee 
OO DCOQOO0OO00O 
+++ Hee ttt 
LU LOL LU Loy Li Ld 
OR- OBDOOIRAIAD 
OOrmrOOAACMO 
OMNAATO SMe 
OO Oerre— rer re~ © 
Om QT UO OO 


oeeé@en32ee0f8f @ @ @ @ 
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OO tee OOF OOO 
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OOO CVOOO0O 
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MUI UU LU bi 
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INDOOR ODA AAD 
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BROWNE DIGIT 1 


TIME OUTPUT VELOCITY ERROR 


e,0000E+00 ~99OUS3ZEt+00 -.24432E-01 © S4750E-02 


NORMAL STOP AT FINAL TIME 
GRAPH TITLED. . BROWNE DIGIT 12 


RUN 1 OUTPUT VS. TIME 
GRAPH TITLED. . BROWNE DIGIT i 
RUN 1 VELOCITY VS.TIME 


GRAPH TITLED. . ge DIGIT 1 


BRC 
RUN EDOQT VS. ERROR 


THE ONE RUN CALLED FOR HAS BEEN COMPLETED. 


STOP 
TIME, 11 MINUTES ANC 14& SECONDS 
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APPENDIX III 


CONTINUOUS PLANT 


The continuous plant which was to be controlled in a sampled- 
data manner consisted of two operational a.c. amplifiers, a d.c. power 
amplifier, and a d.c. torque motor manufactured by the Inland Motor 


Corporation of Pearl River, New York. 






Gee 
Amplifier 






Motor ee 





a.c. 
Amplifiers 


Figure 35 


As is seen in the figure, the control signal to the plant was 
received from the zero-order hold and fed into the first operational 
amplifier. The second operational amplifier transmitted the signal to 
the d.c. power amplifier where the signal size was increased further. 
From the d.c. power amplifier the control signal was fed to the armature 
of the motor. 

The Inland motor has a permanent magnet d.c. field and is armature 
controlled. An inertia disc was attached to the shaft to slow down the 
response of the system. Also on the motor shaft were a tachometer for 
velocity pick-off and a potentiometer which was used for position feed- 


back. The motor was rated at 0.47 ampere and 90 volts. [5 | 
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To provide the needed power to the motor an amplifier whose 
capabilities were greater than the common commercial type was needed. 
The d.c. amplifier used was the same one that had been constructed for 
use on a previous thesis using the same motor. [5 | As seen in Figure 
35 the input to the d.c. amplifier was double ended. The amplifier 
operation is push-pull which reduces third harmonic distortion and 
provides a balanced output. (See Figure 36 for a circuit diagram of 
the d.c. power amplifier.) 

To insure that the d.c. power amplifier was operating correctly 
two tests were performed. The first test ’checked the frequency response 
characteristics of the a.c. and d.c. combination. The second test 
plotted the output vs. the input of the d.c. amplifier alone thus 
checking the saturation characteristic of the amplifier. Test data 
and curves are shown on pages 82-8. It should be mentioned that some 
difficulty was encountered with the d.c. amplifier. The amplifier needs 
a supply voltage of 300 volts with a corresponding current of one 
ampere. Since this power was not available directly in the room where 
the equipment was set up, it was decided to place four portable power 
supplies in parallel. This arrangement proved fairly satisfactory with 
the exception that the power supplied was not regulated. This caused 
difficulty in balancing the amplifiers and once they were balanced they 
did not remain so for any length of time. This balancing problem was 
bothersome but in no way invalidated the results or conclusions of 
this report. 


To determine the transfer function of the amplifier motor 


80 ' 














combination a standard frequency response test was conducted on the 
closed loop system. The input sinusoid was supplied by a Hewlet 
Packard Low-Frequency Generator Model 202A. Unity feedback was em- 
ployed and a Mark II Brush Recorder was used to measure the input and 


output. (See Figure 3/7 ) 














Operational Cae 
aac. Power Motor 
Amplifiers Amplifier 
G(S) 
record record 
Figure 3/ 


The recorded closed loop data was first plotted as a Bode 
diagram. (Magnitude and Phase) The closed loop data was then con- 
verted to an open loop transfer function by using a Nichols Chart and 
Mee cond Bode plot. The open loop transfer function was determined to 

A 


be: G(S) = Scarneo)) with A set at 40. (See pages 85-89 for test 


data, Bode diagrams, and Nichols Chart.) 
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Frequency Response of Amplifiers in Cascade 
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Saturation Characteristic of d.c. Amplifier 
Brequency= 0.7 cps; qs =» 4.39 


input output 
(volts)! (volts) 
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APPENDIX IV 


MEASUREMENT OF QUANTIZATION ERRORS 


This appendix will cover the determination of quantization errors 
in the sampled-data control system. It will only deal with the real 
plant since the main thesis effort was directed to the control of this 
plant. The theory used was described by John B. Slaughter of the Navy 


Electronics Laboratory in the January 1964 IEEE Proceedings on Automatic 
Sontrol . | 6 | 

Quantization, or round off, errors occur because of the inevitable 
conversion operations required when analog and digital devices are 
connected together in a sampled-data closed loop system. Quantization 
is the process of converting a signal in analog form to its digital 


approximation. This is best illustrated by Figure 44. 


Q [x(t] 












ly) 2h 2 





xs) 


Figure 44 


The dashed line in the figure is the desired linear response, while 
the staircase represents the actual output of the quantizer. The 
range of input magnitudes is divided into disjoint intervals (h, ) which 


are not necessarily equal. All magnitudes falling within an interval 


90 





are equated to a single value within the interval. This results in 

the digital approximation to the analog input. It can be seen that the 
maximum round off error of the quantization operation is h/2. (h is 
called the quantizing level.) 

Quantization is non-linear in the sense that a discrete set of 
amplitude levels is produced for some continuous range of input. Be- 
cause of this an analytical expression for the effects of this error 
is difficult to obtain. However, by assuming that the quantizing 
levels are kept small the concepts of probability may be used. The 
quantizer is treated as a summing point with uniformly distributed 
noise introduced at that point. 

In continuing, we will refer to the real plant which was used in 


the sampled-data control problem. 





Figure 45 


In Figure 45 we have the system that was used in this study. The 
sampling period is 0.3 seconds, and the digital controller is designed 


to give a "deadbeat response" to a step input. From z-transform theory 


we have: 
-1 -2 
G(Z) = .368Z +.2642 
-l -2 (1) 
1 - 1.3682 +.368Z 
-1 
2.382 = 2.5622 
D(Z) = eo (2) 
1 + .418Z 
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The equations describing the system may also be written as: 
x(k + 1) = Ax(k) + Dr(k) <3) 


c(k) = Bx(k) (4) 


Now if we decompose our system into its canonical form equations 


2 and 3 may be written by inspection. 








Figure 46 


Canonical form of sampled-data control system 


Thus, if we define our states as shown in Figure 46 we arrive at: 


0 1 0 
n=) Weseesd 833 Beye 
(5) 
SO GE an 
0 
De 42 (6) R= | 1.01 ees 0] (7) 
1 


Before the input, r(k), or the output, c(k), are introduced into 
the digital controller they pass through the A/D converter and are 


changed from analog signals to digital signals. Also, since the 


o2 





digital computer must operate with finite word lengths, round off occurs 
for each word the computer processes or outputs. Therefore quantization 
operations occur and should be accounted for wherever these processes 
happen. This may be done by introducing quartizers in the canonical 


form of the system as shown in Figure 47. 


42 Q 


“a> >> gai *3)- 1711-41 Ore i 
a Q 1.407 





Figure 47 


We now introduce a new vector R(k) which is the difference between 
the quantized and un-quantized terms. The quantized system may now be 


described by: 
x (Ktl)= Ax (kK) + R(k) + Dr(k) (8) 


where the subscript q indicates the quantized system. It may be 
seen by inspecting the above equation that if the un-quantized system 
is stable then the quantized system is also stable. 

To obtain the difference between the quantized and un-quantized 


systems we define: 


e(k) = Eel - x(k) (9) 
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from which: 


ei gene CK) en) (10) 

cg(k) - e(k) = Be(k) (11) 
Thus in our problem: 

e, (k+l) = e, (k) (12) 


e, (k+1)= ~-83le, (k) + .833e, (k) - - 349e,, (k) + R, (k) (13) 
e, (kt+1)= ~1,0le, (k) ~ 1.365e, (k) - -424e, (k)+R, (Kk) (14) 


We must now find an upper bound on R and hence one. From this we 
may obtain the maximum output error due to quantization. Since R. (k) 
is the difference between the state vector component, x, (k+l), obtained 
in the quantized and un-quantized systems, an upper bound on R. (k) is 
the greatest error due to quantizing in one sampling period. 


In our system we have: 


R, (k) = 0 (15) 
R,(k) £ | 42h, + .2Ih, + h,h,| 2 H (16) 
2 Deere? een eagle 2 
< ca 
R,(k) € Ih + <5h, H, (17) 





where h,/2 is the A/D quantizing error, h,/2 is the computer 
round off error, and ha/2 is the D/A quantizing error 
To find an upper bound on e(k) we must convert equations 12, 13, and 
14 in their z-transform equivalents. 


Thus: 


252) = ze, (z) ° (18) 





all 


oa a si pe 
e,(z)= -.83lz “e,(z) + .833z “e,(z) - .349z e,(z) + -—— (19) 
2 1 2 3 int 

zi a Es Haz 
e.(z)= -1.01z e.(z) - 1.3652 e.@) = 74242 e.(z)+ —— (20) 
3 1 2 5 1-27 


Referring to equations 7 and 11 it is seen that to find the overall 
system error due to quantization we only have to solve for e, (2) and e, (Z). 


Also, the steady state error e, equals the steady state error e That 


1 2° 


ms< lim Gee ye (z) = lim Gee oN (z) 
iL 2 
Z>l1 z>l 


Therefore after solving equations 18, 19 and 20 simultaneously and 


applying the final value theorem we arrive at: 


1.4248, - 3491, 


e1 (co) = ey (oO) = (oa (21) 
195 
and since in the hardware used hy = 2.44 millivolts, 
h, =), 024¢. mr Mlavotes sand hy = 2.44 millivolts 
then 
er (oO) = en (oo) = 15.3 millivolts (22) 


Therefore Sa) - c(k) = 36.35 millivolts 

This difference between the quantized and unquantized systems was less 
than 1% of the final value of the system for a five volt step which was 
ordinarily used. Therefore no problems were encountered with quantization 


errors in the control of the real system. 


oh, 





APPENDIX V 


PROGRAM FLOWCHART SYMBOLS 


denial 


Represents 


Processing 


External Function 


Terminal or Connector 


Decision 


Indexing Operation 


Predefined Process 


Jump 


96 





APPENDIX VI 


PROGRAMS FOR THE CDC-160 COMPUTER 


PROGRAM 1 
Location Contents Comments 
1000 7500 Call A/D 
1001 1401 
1002 7600 Read into A 
1003 L000 Store in 0000 
1004 7500 Call D/A 
1005 2h01 
1006 7304 Output (0000) 
1007 0001 
1010 6510 Jump to 1000 
1011 6411 
1012 0000 


oy 


eS «& 
—- S 
= > 
> 
* aa _ 





Location 


1000 
1001 
1002 
1003 
1004 
1005 
1006 


1007 


1010 


1011 


1012 


1013 
1014 
1015 
1016 


1017 


1020 
1021 
1022 
1023 
1024 


1025 


Contents 


7500 
1401 
7600 
Loé2 
3460 
LO65 
3063 
6315 


2065 
3463 
6211 
2062 
4060 
7500 
2401 
7306 


0061 
6521 
éhe2 
7071 
7072 
0060 


Program 2 


98 


Comments 


Call A/D 


Read into A 


Store 'New value’ 
in 0062 

Subtract ‘Old 
Value' 

Store 'Diff' in 
0065 


Add Increment 

If sum negative, 

jump to 1024 

Load 'Diff' into A 
Subtract Increment 

If difference positive, 
jump to 1023 

Load 'New Value! 

Store as 'Old Value! 


in 0060 
Call D/A 


Output (0060) 


Jump to 1000 


Jump to Negative Step 
Solution (0071) 
Jump to Positive Step 
Solution (0072) 


PROGRAM 2 (CONTINUED) 


LOW CORE ALLOCATION 


Location Contents Comments 

0060 — 'Old Value' 
0062 ---- 'New Value' 
0063 0010 Increment 

0065 —e 'Diff' 

0071 T7100 Dummy start of 


0072 


T(00 


2, 


Negative Step 
Solution 


Dummy start of 
Positive Step 
Solution 








Location 


1030 
1031 
1032 
1033 
1034 
1035 
1036 
1037 


1040 
1041 
1042 
1043 
1044. 
1045 
1046 
1047 


1050 
1051 


1052 


PROGRAM 3 


Contents 


7500 
2401 
TALT 
0021 
2050 
3450 
4050 
5450 


2050 
Bu51 
6002 
6504 
1900 
2401 
7304 
0023 


T{00 
0020 


Comments 


Call D/A 


Output (0020) 


zeroize counter 


Add one to 
counter 


Subtract counter 
total from 
preset total 


If difference = QO, 
jump to 1lokh 

If difference not 
= O, jump to 1057 


Call D/A 


Output (0022) 


Halt 





. 








PROGRAM 3 (CONTINUED ) 


LOW CORE ALLOCATION 


Location Contents Comments 
0020 0631 Converted as 
-k volts 
0022 7146 Converted as 
-6 volts 
0050 seme Counter 
0051 3403 Preset total for 


0.1 second delay 


Tor 








Location 


1000 
1001 
1002 
1003 
1004 
1005 
1006 


1007 


1010 
1011 
1012 
1013 
1014 
1015 
1016 


1017 


1020 
1021 
1022 
1023 
1024 


1025 


PROGRAM 4 


Contents 


1300 
1401 
7600 
4062 
3460 
4065 
3063 
6315 


2065 
3463 
6211 
2062 
LO60 
7500 
2ho1 
7306 


0061 
6521 
6422 
7071 
1072 
0060 


102 


Comments 


Sense Loop 


(See PROGRAM 2) 





Location 


1026 


1027 


1030 
1031 
1032 
1033 
1034 
1035 
1036 
1047 


1040 
1042 
1042 
1043 
1044 
1045 
1046 
1047 


1050 
1051 


1052 


PROGRAM 4 (CONTINUED ) 


Contents 


2060 
3014 


4020 
2060 
3015 
Lo2e 
71500 
ehol 
(627, 
0021 


2052 
3452 
Lo52 
2050 
3450 
4O50 
5450 


2050 


5451 
6002 
6504 


103 


Camments 


NEGATIVE STEP SOLUTION 
Load ‘Old Value' 


Add H° 


Store adjusted H° 
in 0020 

Load ‘Old Value! 

Add H} 

Store adjusted H? 


in 0022 
Call D/A 


Output (0020) 


Start Delay 


(See PROGRAM 3) 








Location 


1053 
1054 
1055 
1056 
1057 


1060 
1061 
1062 
1063 
1064 
1065 
1066 
1067 


1070 
1071 
1072 
1073 
1074 
1075 
1076 
LOTT 


PROGRAM 4 (CONTINUED) 


Contents 


5452 
2052 
3453 
6002 
6514 


7500 
2401 
7304 
0023 
1073 
0020 
0022 
2060 


3016 
Lo2k 
2060 
3017 
4026 
7500 
2ho01 
TOEN 


04 


Comments 


Call D/A 


Output (0022) 


Jump to 1130 


POSITIVE STEP SOLUTION 
Load ‘Old Value! 


Ada H° 

Store adjusted He 
in 002k 

Load ‘Old Value’ 
Ada H+ 

Store adjusted H* 


in 0026 
Call D/A 


Output (0024) 





Location 


1100 
pbkeyl 
1102 
1103 
1104 
1105 
1106 


1107 


1110 
abut 
IL 
aR ils, 
1114 
1115 
1116 


1120 
1121 


iv2e2 


12h 


PROGRAM 4 (CONTINUED) 


Contents 


0025 
2052 
3452 
hos52 
2050 
3450 
4O50 
5450 


2050 


3u51 


6504 
5452 
2052 
3453 
6002 


6514 
1900 
2hbo1 
7304 
0027 


LU5 


Comments 


Start Delay 


(See PROGRAM 3) 


Call D/A 


Output (0026) 





Peery ee 


. 





Location 


125 
1126 


ll27 


1130 
piow 
1132 
1153 
1134 
1135 
1136 
1137 


1140 
1141 
1142 
1143 


1145 
1146 


1147 


1150 


PROGRAM & (CONTINUED ) 


Contents 


7074 
0024 


0026 


2052 
3452 
4052 
2050 
3450 
Lo50 
5450 


2050 


6002 
6504 
5452 
2052 
3453 
6002 
6514, 


1900 


106 


Comments 


Jump to 1130 


Start Delay 


(See PROGRAM 3) 


BOOTS TRAP 





tI 





Location 


1151 
1152 


dL a 


PROGRAM 4 (CONTINUED ) 


Contents 


1401 
7600 
4060 


1154 


0014 
0015 
0016 
0017 
0020 
0022 
002k 
0026 
0051 


0053 


0060 
0062 
0063 
0065 
0070 
0071 


0072 


0073-0074 


1070 


LOW CORE ALLOCATION 


0005 


1000 
1026 
1067 
1130 


107 


Comments 


Read into A 
Store as ‘Old 


Value! 
Jump to 1000 


H° (-) 

H* (-) 

H° (+) 

H? (+) 

Adjusted H®° (-) 
Adjusted H? (-) 
Adjusted H° (+) 
Adjusted H? (+) 
Preset total for 
O.1 second delay 
Number of passes 
through 0.1] second 
delay 

'Old Value' 

'New Value! 


Increment 


‘Dis 





Location 


1000 
1001 
1002 
1003 
1004 
1005 
1006 


1007 


1010 
1011 
1012 
1013 
1014 
1015 
1016 


1017 


1020 
1021 
1022 
1023 
1024 


1025 


PROGRAM 5 


Contents 


1500 
1401 
7600 
4062 
3460 
LO65 
3063 
6315 


2065 
3463 
6211 
2062 
4O60 
7500 
201 
7306 


0061 
6521 
6420 
7011 
7012 
0060 


108 


Comments 


Sense Loop 


Jump to Negative 
Step Solution 
Jump to Positive 
Step Solution 








Location 


1026 


1027 


1030 
1032 
1032 
1053 
1034 
1055 
1036 
1037 


1040 
1042 
1042 
1043 
1044 
1045 
1046 
1047 


1050 


1051 


PROGRAM 5 (CONTINUED ) 


Contents Comments 


1500 
1401 


7600 
3460 
Look 
0101 
7000 
0031 
0020 
002k 


0050 
0101 
7000 
0041 
0030 
0002 
0030 
2030 


3060 
4030 


Negative Step 
Solution 
Call A/D 


Read into A 
Subtract 'Old 
Value! 

Store E° 


Multiply E° 
by 100 K® 


Divide product 
by 100 


Load H°? 


Add 'Old Value’ 


Store adjusted H° 


109 





Location 


1052 
1053 
1054 
1099 
1056 
1057 


1060 
1061 
1062 
1063 
1064 
1065 
1066 
1067 


1070 
1071 
1072 
1073 
1074 


PROGRAM 5 (CONTINUED ) 


Contents 


1900 
2ho1 
1393 
0031 
2052 
3452 


4o52 
2050 
3450 
4050 
5450 
2050 
BU51 
6002 


6504 
5452 
2052 
3453 
6002 


110 


Comments 


Call D/A 
Output adjusted 
He 


Delay 





»>esstec ee 









Location 


1075 
1076 
1077 


1100 
1101 
1102 
1103 
1104 
LLO5 
1106 


1107 


1110 
1111 
1112 
1113 
1114 
1115 
1116 
a ey 


1120 


PROGRAM 5 (CONTINUED ) 


Contents 


6514 
7500 
1401 


7600 
3460 
Look 
0101 
7000 
0031 
0022 


0024 


0030 
0101 
7000 
0041 
0030 
0002 
0050 
2030 


3060 


Comments 


Call A/D 


Read into A 
Subtract 'Old 
Value' 

Store Et 


Multiply E? 
by 100K2 


Divide product 
by 100 


Load H* 


Add ‘Old Value' 





Location 


Bei. 
1122 
1123 
1124 
1125 
1126 


1127 


1150 
1131 
1132 
1155 
1134 
1135 
1136 
1137 


1140 
1141 
11he 
1143 
Whb 


PROGRAM 5 (CONTINUED ) 


Contents 


4030 
1900 
2h01 
1305 
0031 
O13 
00350 


1500 
1401 
7600 
3460 
Loak 
ahah 
Lo2k 


0101 


7000 
0051 
0020 
0024 


0030 


Mie 


Comments 


Store adjusted H? 


Call D/A 


Output adjusted Ht 


Jump to second 
delay 


Positive Step 
Solution 
Call A/D 


Read into A 
Subtract ‘Old 
Value' 

Store E° 


Complement E° 


Multiply complemented 
E° by 100 kK®° 





Location 


1145 
1146 
1147 


1150 
1151 
1152 
1153 
1154 
1155 
1156 
L157 


1160 
1161 
1162 
1163 
1164 
1165 
1166 
1167 


1170 


PROGRAM 5 (CONTINUED) 


Contents 


0101 
[000 


oo41 


0030 
0002 
0030 
24.30 
4030 
2030 
3060 
4030 


1300 
2401 
1397 
0031 
2052 
3452 
Lo52 


2050 


3450 


113 


Comments 


Divide product 
by 100 


Complement 
dividend 

Store H° 

Load H° 

Add ‘Old Value’ 


Store adjusted H° 


Call D/A 


Output adjusted H° 


Delay 





Location 


1171 
1172 
1173 
1174 
1175 
1176 
1177 


1200 
1201 
1202 
1203 
1204 
1205 
1206 


1207 


1211 
a =a 2 
gb i 
1214 


1215 


PROGRAM 5 (CONTINUED) 


Contents 


4O50 
5450 
2050 
Bh51 
6002 
6504 
5452 


2502 
3453 
6002 
6514 
7500 
1401 
7600 
3460 


hook 
ehok 
Look 
0101 
7000 


0031 


114 


Comments 


Call A/D 


Read into A 
Subtract ‘Old 
Value! 

Store Et 


Complement E21 


Multiply complemented 
E2 py 100K2 


Location 


1216 


Lea 


1220 
1221 
1222 
1223 
1224 
1225 
1226 


1227 


1230 
1231 
1232 
1253 
1234 
1235 
1236 
1257 


12h0 


PROGRAM 5 (CONTINUED) 


Contents 


0022 


002k 


0030 
0101 
7000 
ook 
0030 
0002 
0030 
2430 


4030 
2030 
3060 
4030 
71500 
2ho1 
1503 
0031 


(013 


5 


Comments 


Divide product 
by 100 


Complement 
dividend 


Load H? 
Add 'Old Value! 
Store adjusted H? 


Call D/A 


Output adjusted H? 


Jump to second 


delay 





Location 


124) 
1242 
1243 
12h 
12k5 
1246 
1247 


1250 
1251 
1252 
1253 
1254 
1255 
1256 
Je all 


1260 
1261 
1262 
1263 
1264, 
1265 
1266 


PROGRAM 5 (CONTINUED) 


Contents 


0030 
2052 
3452 
L052 
2050 
3450 
LO50 


5450 
2050 
3451 
6002 
6504. 
5452 
2052 


3453 


6002 
6514 


7500 
1401 


Lo60 


7010 


116 


Comments 


Second Delay 


Bootstrap 


PROGRAM 5 (CONTINUED ) 


LOW CORE ALLOCATION 


Location Contents Comments 

0000 wecn Address of 
Subroutine ARITH 

0002 O144 

0010 1000 Address of start 
of program 

0011 1026 Address of negative 
step solution 

0012 1130 Address of positive 
step solution 

0013 1164 Address of second 
delay 

0020-21 ---- K° 

0022-23 ---- K* 

0024-25 ---- Current error 

0030-31 sne- Current H 

0051 3403 Preset total for 
O.1 second delay 

0053 0005 Number of passes 
through 0.1 second 
delay 

0060 oe 'Old Value’ 

0062 on-- 'New Value’ 

0063 ---- Increment 

0065 ---- 'Dift' 





PROGRAM 6 


Due to the basic similarity of this program to the 
preceding programs, the comments presented here will be 
general in nature. 


Location Contents Comments 

1000 7500 SENSE LOOP 
1001 1401 

1002 7600 

1003 Lo6h 

1004 7500 

1005 1402 

1006 7600 

1007 LLo6é5 

1010 2064 Sense error in 
1011 3461 System 1 

1012 4067 

1013 3063 

1014 6351 If error, set 
1015 2067 BE} Flag 

1016 3463 

1017 6250 

1020 2064 

1021 4o61 

1022 6101 

1023 4066 

1024 7500 If no error, D/A 
1025 2ko1 System 1 

1026 1353 

1027 0067 

1030 2065 Sense error in 
1031 3462 System 2 

1032 4067 

1033 3063 

1034 6337 If error, set 
1035 2067 E° Flag and 
1036 3463 jump to solution 
1037 6236 routine 

1040 2065 

1041 4o62 

1042 6101 

1043 ho66 

1044 7500 


118 





Location 


1045 
1046 
1047 


1050 
1051 
1052 
1053 
1054 
1055 
1056 
1057 


1060 
1061 
1062 
1063 
1064 
1065 
1066 
1067 


1070 
1071 
1072 
1073 
1074 
1075 
1076 
1O77 


1100 
1101 
1102 
1103 
1204 
1105 
1106 
1107 


1110 
Lash) 
1lie 
1113 


PROGRAM 6 (CONTINUED ) 


Contents 


2ho2 


liebe, 
0067 


2040 
6123 
6122 
2050 
3450 
4050 
5450 
2050 


3451 
6002 
650h4 
6463 
6564 
SLO 
6536 
20ho0 


3401 
Lolo 
6542 
Shh. 
7010 
20h1 
3401 
Lok 


7010 
0066 
0000 
0000 
0000 
2067 
3467 
4067 


2040 
6003 
6304 
6205 


pee, 


Comments 


If no error, 
D/A System 2 


If E* Flag set, 
jump to Solution 
Routine 


Basic Delay 


Return to start 
of Sense Loop 


SOLUTION ROUTINE 
E? Flag test 


E* Flag test 











Location 


1114 
1115 
1116 
1117 


1120 
1121 
1122 
1123 
1124 
1125 
1126 
1127 


1130 
1131 
1132 
1133 
1134 
1135 
1136 
1137 


1140 
1241 
1242 
1143 
1244 
1145 
1146 
1147 


1150 
1151 
1152 
1153 
1154 
pp 
1156 
57 


1160 
1161 
1162 


PROGRAM 6 (CONTINUED) 


Contents 


7101 
13570 
7101 
1250 


eoke 
6002 
6137 
7500 
1401 
7600 
3461 
4020 


2420 
4o20 
0101 
7000 
0031 
0020 
0022 
0016 


0101 
7000 
Ook.1 
0016 
0004 
0016 
216 
3061 


4o16 
7500 
2ho1 
7305 
0017 
Si.ke 
7101 
1370 


0016 
ae) 
2044 


120 


Comments 


If = -l1, jump to 
Negative Step 
Solution for System 1 


Positive Step 
Solution for 


System 1 
If J #0, jum to 1161 


Compute and output 
H°(+) for System 1 


J=J+#i1 


Jump to E* Flag test 


T=I+1 





Location 


1163 
1164 
1165 
1166 
1167 


1170 
1171 
1172 
1173 
1174 
1175 
1176 
1177 


1200 
1201 
1202 
1203 
1204 
1205 
1206 
1207 


1210 
Lei) 
1212 
Lae 
1e14 
1215 
1216 
Ley 


1220 
leer. 
1222 
Tae5 
1224 
1225 
1226 
1227 


1230 
1231 


PROGRAM 6 (CONTINUED) 


Contents 


3403 
6506 
20ke 
3401 
6003 


7101 
1760 
7500 
1401 
7600 
3461 
4066 
eh66 


LO66 
3420 
6303 
7101 
1205 
0101 
7000 
0031 


0066 
0024 
0016 
0101 
7000 
ook1 
0016 
ooo4 


0016 
2416 
3061 
4016 


7500 
eho1 


1510 
0017 


Sho 
2043 


zal 


Comments 


If I-# (0003), 
jump to E* Flag test 


If J #1, jump to 
E1R Flag test 


Compute and output 
H?(+) for System 1 


J=J +1 





Location 


1232 
1233 
1234 
1235 
1236 
1237 


12h0 
12h1 
12he 
123 
12h 
125 
126 
1247 


1250 
1251 
1252 
1253 
1254 
1255 
1256 
1257 


1260 
1261 
1262 
1263 
1264 
1265 
1266 
1267 


1270 
1271 
1272 
vero 
1274 
1275 
1276 
le77 


PROGRAM 6 (CONTINUED ) 


Contents 


343 
Los 
7101 
1370 
0016 
0000 


27 


Jump to E* Flag test 


Negative Step 
Solution for 
System 1 


(This solution is 
handled exactly as 
is its positive 
counterpart. Thus, 
no comments will be 
made. ) 





Location 


1300 
1301 
1302 
1305 
1304 
1309 
1306 
1307 


1310 
1311 
1312 
1313 
1314 
1315 
1316 
1317 


1320 
1321 
1322 
1323 
1324 
1325 
1326 
1327 


13350 
1331 
13532 
L559 


1359 
1336 
1337 


1340 
1341 


1343 
1344 
1345 
1346 
1347 


PROGRAM 6 (CONTINUED) 


Contents 


2ho1 
1309 
0017 
54h2 
7101 
1370 
0016 
Shh 


2043 


3403 
6506 
2042 
3401 
6003 
7101 
1760 


7500 
1401 


4066 
3420 
6303 
7101 
1530 


0101 
7000 
0031 
0066 
002k 
0016 
0101 
7000 


ook] 
0016 
ooo4 
0016 
2016 
3061 
4o16 
7500 


iz 





Location 


1350 
1351 
Lo 
1393 
1354 
1355 
1356 
1397 


1360 
1361 
1362 
1363 


1365 
1366 
1367 


1370 
1371 
1372 
Bi fe) 
1374 
159 
1376 
alo | 


1400 
1401 
1402 
1403 
hob 
1405 
1406 
1407 


1410 
1411 
1412 
1413 
1414 
1415 
1416 


PROGRAM 6 (CONTINUED ) 


Contents 


eho1 
7310 
0017 
5hke 
2043 
34h 
Loks 
7101 


1370 
0016 
0000 
0000 
0000 
0000 
0000 
0000 


2047 
6003 
6304 
6205 
7101 
1650 
7101 
1525 


20k 
6002 
6137 


7500 
1402 


124 


Comments 


E? Flag test 


If = 0, jump to 


delay 


If = -l, jump to 
Negative Step 
Solution for 
System 2 


Positive Step 

Solution for 

System 2 

If L #4 0, jump to 
hh 


1441 


Compute and output 
H°(+) for System 2 











Location 


1417 


1420 
1427 
14202 
1423 
1heh 
1425 
1426 
1427 


1430 
1431 
1432 
1433 
1434 
1435 
1436 
1437 


rho 
hh 
14he 
1443 
Tybhk 
Ths 
1446 
1447 


1450 
1451 
1452 
1453 
1454, 
1455 
1456 
1457 


14.60 
1461 
1462 
1463 
1464 
1465 


PROGRAM 6 (CONTINUED) 


Contents 


0016 


0101 
7000 
0041 
0016 
0004 
0016 
2416 
3062 


4016 
7500 
aloe 
7305 
0017 
Shh 
7101 
1650 


0016 
5uks 
e20ks 
3402 
6506 
e20k4 
3401 
6003 


7101 
2040 
7500 
1402 
7600 
34.62 
4o66 
466 


3430 
6303 
7101 
1464 
0101 
7000 


LZ 


Comments 


L.= Lat 


Jump to delay 


K=K+1 


If K # (0002), 
jump to delay 


If L #1, jump to 
ER Flag test 


Compute and output 
H1(+) for System 2 


————— 


A 
Jaa 
: 
® 
: 
or 
't4 
, 
‘ 
«@ 
@ 
: 
as < 
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Location 


1466 
1467 


1470 
1471 
1472 
14-73 
1474 
1475 
1476 
1477 


1500 
1501 
1502 
1503 
1504 
1505 
1506 
1507 


1510 
1511 
1512 
1513 
1514 


. ali lp 


1516 
1517 


1520 
1521 
1522 
1523 
1524 
1525 
1526 
1527 


1530 
DSN 
1552 
193) 
1534 


PROGRAM 6 (CONTINUED) 


Contents 


0031 
0066 


0034 
0016 
0101 
7000 
ook1 
0016 
0004 
0016 


2476 
3062 
4o16 
71500 
elioe 
7310 
0017 
Shu 


2045 
3hh5 
Loks 
7101 
1650 
0016 
0000 
0000 


0000 
0000 
0000 
0000 
0000 
204k 
6002 
6135 


7500 
1402 
7600 
3462 
4030 


126 


Comments 
L=L+t+1l1 
K = 0 


Jump to delay 


Negative Step 
Solution for 
System 2 


(As before, no 
comments will be made 
for the negative 
step solution.) 





= i) OO iv 


& = 


= 





Location 


1559 
1536 
L557 


1540 
1541 
1542 
1543 
1544 
1545 
1546 
1547 


1550 
1551 
1552 
155 
1554 
1)? 
1556 
1557 


1560 
1561 
1562 
1563 
1564 
1565 
1566 
1567 


1570 
1571 
1572 
1573 
1574 
15D 
1576 
L5TT 


1600 
1601 
1602 
1603 


PROGRAM 6 (CONTINUED ) 


Contents 


0101 
{000 
0031 


0030 
0032 
0016 
0101 
T7000 
oOk1 
0016 
OOO 


0016 
2016 
3062 
4016 
7500 
2ho2 
7305 
0017 


Sukh 
7101 
1650 
0016 
Suk5 
20k5 
3402 
6506 


20k 
3401 
6003 
7101 
2040 
7500 
1402 
7600 


3462 
4O66 
3430 
6303 


ney 








Location 


1604 
1605 
1606 
1607 


1610 
1611 
1612 
1613 
1614 
1615 
1616 
1617 


1620 
1621 
1622 
1623 
1624 
1625 
1626 
1627 


1630 
1631 
1632 
1633 
1634 
1635 
1636 
1637 


1640 
1641 
1642 
1643 
1644 
1645 
1646 
1647 


1650 
1651 
1652 
1653 


PROGRAM 6 (CONTINUED) 


Contents 


7101 
1606 
0101 
7000 


0031 
0066 
0034 
0016 
0101 
7000 
ook1 
0016 


ooo4. 
0016 
2016 
3062 
4016 
7500 
2hoe 
7310 


0017 
Suuy 
20k5 
3445 
Loks 
7101 
1650 
0016 


0000 
0000 
0000 
0000 
0000 
0000 
0000 
0000 


2050 
3450 
LO50 
5450 


128 


Comments 


Delay 








Location 


1654 
1655 
1656 
1657 


1660 
1661 
1662 
1663 
1664 
1665 
1666 
1667 


1670 
1671 
1672 
1673 
1674 
1675 
1676 
1677 


1700 
1701 
1702 
1703 
1704 
1705 
1706 
1707 


1710 
1711 
1712 
1713 
1714 
1715 
1716 
slvaly4 


1720 
jen 
1722 


PROGRAM 6 (CONTINUED ) 


Contents 


2050 
3451 
6002 
6504 


200 
6133 
7500 
1401 
7600 
Lo6k 
3461 
4067 


3063 
6316 
2067 
3463 
6215 
2064 
4O61 
6101 


ho66 


7500 
201 
7303 
0067 
6107 
0066 
SLO 


6104 
20ho 
3401 
Lolo 
2041 
6133 
7500 
1402 


7600 


4O65 
3462 


T29 


Comments 


If E? Flag #1, 
sense error for 
System 1 


If E* Fleg = 0, 
sense error for 
System 2 











Location 


1723 
172k 
1725 
1726 
1727 


1730 
1731 
1732 
> 
1734 
1755 
1736 
1737 


PROGRAM 6 (CONTINUED) 


Contents 


O67 
3063 
6316 
2067 
3463 


6215 
2065 
4o62 
6101 
4066 
1500 
ehoe 
7303 


0067 
7110 
0066 
Shb7 
7105 
2041 
3401 
Lok1 


7010 
1105 
0000 
0000 
0000 
0000 
0000 
0000 


2046 
6132 
eoke 
Zhe 
hoke 
203 
bie) 
Loks 


20hko 
Zuo 


1 30 


Comments 


Jump to start of 
Solution Routine 


RETURN ROUTINE 

If ER Flag # 0, 

jump to ramp solution 
routine (not included) 


J =0 
I=0 
E? Flag = 0 














Location 


1772 
1TT3 
1774 
17719 
gir 
L7TT 


2000 
2001 
2002 
2003 
2004 


2006 
2007 


2010 
2011 
2012 
2013 
2014 
2015 
2016 
2017 


2020 
2021 
2022 
2024 
202k 
2025 
2026 
2027 


2030 
2031 
2032 
20353 
2034 
2035 
2036 
2037 


20k0 


PROGRAM 6 (CONTINUED) 


Contents 


holo 
20h] 
6112 
7500 
1401 
7600 


O61 
7500 
1ho02 
7600 
4062 
7011 


7500 
1401 


7600 
ho61 
7010 
7101 
1105 
0000 
0000 
0000 


Comments 


If E* Flag #0, 
bootstrap System 1 
and jump to Solution 
Routine 


Bootstrap Systems 
1 and 2, and jump 
to Sense Loop 








Location 


2041 
20he2 
2043 
20k 
2045 
2046 
20K7 


2050 
2051 
2052 
2095 
2054 
2059 
2056 
2057 


2060 
2061 
2062 
2063 
2064 
2065 
2066 
2067 


2070 
2071 
2072 
2075 
2074 


0000 


0001 
0002 
0003 
0004 
0010 
0011 


PROGRAM 6 (CONTINUED) 


Contents 


6132 
e0kk 


7600 
4062 
7011 
7500 
1402 


7600 
4062 
7010 
7101 
4000 


LOW CORE ALLOCATION 


0001 


O14 
1105 
1000 


12 


Comments 


If E“R Flag # 0, 


jump to ramp 
solution routine 


=O 
K = 0 
E- Flag = 0 


If E* Flag # 0, 
bootstrap System 2 
and jump to 
Solution Routine 


Bootstrap Systems 
to Sense Loop 


Address of Sub- 
routine ARITH 


Q 
Vig 





Location 


0016-17 
0020-21 
0022-23 
0024-25 
0030-31 
0032-33 
0034-35 
ooo 
Ook1 
ook2 
0043 
ooky 
OO45 
0046 
0047 


0050 
0051 


PROGRAM 6 (CONTINUED) 


Contents 


0675 


0100 


Comments 


H's for output 
E(1) 
100K° (1) 
100K?(1) 
E(2) 
100K° (2) 
100K (2) 
E? Flag 
E* Flag 
J 

I 

L 


K 
E7R Flag 
E°R Flag 


Delay counter 
Preset total 
for 25 milli- 
second delay 


Old Value (1) 
Old Value (2) 
Increment 

New Value (1) 
New Value (2) 
Temp Storage 


Note: Cells 0070 through 0077 and cell 0007 are reserved for 


Subroutine ARITH. 


Leys 


bn) 
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